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Captures of a particles on the proton-richest Barium isotope, ^"^"Ba, have been studied in order 
to provide cross section data for the modelling of the astrophysical 7 process. The cross sections of 
the ^^''Ba(a,7)'^^*Ce and ^''°Ba(a,n)^''^Ce reactions have been measured with the activation tech- 
nique in the center-of mass energy range between 11.6 and 16MeV, close above the astrophysically 
relevant energies. As a side result, the cross section of the ^■^^Ba(a,n)^'^^Ce reaction has also been 
measured. The results are compared with the prediction of statistical model calculations, using 
different input parameters such as a-|-nucleus optical potentials. It is found that the (a,n) data 
can be reproduced employing the standard a+nucleus optical potential widely used in astrophysical 
applications. Assuming its validity also in the astrophysically relevant energy window, we present 
new stellar reaction rates for ^^''Ba(a,7)^^''Ce and ^■^^Ba(a,7)'^^^Ce and their inverse reactions cal- 
culated with the SMARAGD statistical model code. The highly increased ^'^®Ce(7,Q)^'^^Ba rate 
implies that the p nucleus ^^'^Ba cannot directly receive contributions from the Ce isotopic chain. 
Further measurements are required to better constrain this result. 

PACS numbers: 26.30.-k,25.55.-e,27.60.-(-j 



I. INTRODUCTION 

The synthesis of chemical elements heavier than iron 
is a complex problem in nuclear astrophysics [H, |3] • The 
origin of about 99% of these heavy elements observed 
in the Solar System can be rather well explained by the 
classical slow and rapid [il] neutron capture processes 
although details of these processes (especially the r pro- 
cess) are still not completely known. The missing 1 % 
comprises those heavy, proton rich isotopes which can- 
not be created in neutron capture reactions. These are 
the so-called p isotopes and in general their hypothetical 
production mechanism is referred to as the astrophysical 
p process @. 

Since no single process has been found so far which 
could explain the observed abundances of all the p iso- 
topes, the p process may actually involve several sub- 
processes, each of which may have its contribution to the 
production of at least a sub-group of p isotopes. The 
possibilities include, e.g., the rp process @, the v pro- 
cess [7], and the recently studied lyp process The 
most important process is perhaps the so-called 7 pro- 
cess, which - as its name suggests - proceeds through 
7-induced reactions 

If heavy seed isotopes produced earlier by the s or r 
processes are available in a hot stellar environment, then 
the high energy thermal photons may induce (7,n) re- 
actions driving the material towards the proton-rich iso- 
topes. For more and more neutron deficient nuclei, also 
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charged particles can be emitted. Then the (7,p) and 
(7, a) reactions compete with the (7,n) reactions and in- 
fluence the reaction path. The required high temper- 
atures (of the order of a few GK) are encountered in 
evolved massive stars either during the supernova explo- 
sion or shortly before that. Therefore, the 7 process is 
thought to take place in core-collapse supernova explo- 
sions, although other possible sites have also been sug- 
gested [13. 

When the whole mass region of p isotopes is considered, 
the modeling of the 7 process requires a huge reaction 
network of tens of thousands of reactions. The stellar 
rates of these reactions are required input to the network 
to calculate the p isotope abundances. So far, no 7 pro- 
cess calculation in massive stars reproduced well the ob- 
served abundances of the p isotopes simultaneously over 
the whole mass region. Since the reaction rates are typ- 
ically taken from largely untested theoretical cross sec- 
tions, the failure of the models may at least in part be 
attributed to incorrect reaction rates. Therefore, the ex- 
perimental veriflcation of the calculated cross sections is 
of high importance to provide better input for the net- 
work calculations and increase their reliability. 

In order to reduce the complications caused by the 
thermal population of target states in the astrophysi- 
cal plasma, it is preferable to study experimentally the 
inverse capture reactions instead of the 7-induced reac- 
tions [ill, • Comparison of the available experimen- 
tal data with the theoretical cross sections (which are 
typically obtained from Hauser-Feshbach type statistical 
model calculations) shows that for (p,7) and (n,7) reac- 
tions the models are able to reproduce the measured cross 
sections within about a factor of two. This can be con- 
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FIG. 1. (color online) The investigated reactions and the decay of the reaction products. The three studied reactions are 
indicated by the thick grey arrows while the decays which have been used for the analysis are shown by the blue arrows. All 
the decay data are taken from refs. [l^ - [20l ) with the exception of the ^^^Ce"* half-life, where the newly determined value is 
used fH. 



sidered as the general precision of the statistical models 
for these reactions, some fine-tuning of the parameters, 
however, being stiU possible [l^l- For (a, 7) reactions, 
on the other hand, much larger deviations can be found. 
Moreover, there are only a handful of isotopes for which 
{a, 7) cross section data are available in the mass range 
relevant for the 7 process. This severely limits the possi- 
bility of testing the model calculations (for a recent list 
of measured reaction see [l^)- Therefore, measurements 
of a capture reaction cross sections in the relevant mass 
and energy range may help increase the reliability of sta- 
tistical model calculations and lead to a better modeling 
of the 7 process. 

Not all the reactions in a 7 process network have equal 
importance for the resulting p isotope abundances. By 
carrying out a 7 process network calculation and vary- 
ing the involved reaction rates, one can pinpoint those 
reactions to which the calculations show the highest sen- 
sitivity. Such a study has been carried out recently and 
a list of isotopes has been suggested on which proton or 
a capture cross sections should be measured |15| . From 
the 15 naturally occurring isotopes for which high pri- 
ority has been assigned to (a, 7) measurements, the p 
isotope of the element barium, ^'^"Ba, has been selected 
for the present study. No experimental data is available 
yet in the literature for a-induced reactions on ^'^"Ba. 

It can be shown (e.g., [l^ and also the discussion in 
Sec. fVl below) that the (a, 7) cross sections at astrophys- 
ical energies almost exclusively depend only on the a 
widths. Therefore they show a strong sensitivity to the 
a-l-nucleus optical potential in statistical model calcu- 
lations. Global optical potentials designed for a wide 
range of isotopes are typically used in astrophysical ap- 
plications. Several global potentials are available in the 
literature and different potentials lead to largely differ- 
ent cross sections. The comparison with experimental 



data can aid the selection of the best optical potential 
for the description of low energy a capture reactions in 
the heavy element range and thus the collection of more 
experimental data at low energy is useful also for a better 
understanding of the optical potential. 

In the following, we present further details of the inves- 
tigated reactions in Sec. |lll The experimental procedure 
is thoroughly explained in Sec. lIIII and the resulting cross 
sections are presented in Sec. IIVI The impact of the data 
on predictions of astrophysical reaction rates for the reac- 
tions i30Ba(a,7)i34Ce and ^^'^Ce{j,ay^°Ba. is discussed 
in Sec. |Vl Tables of the recommended reaction rates and 
their fits are presented in Sec. I VII along with some con- 
cluding remarks. 



II. THE INVESTIGATED REACTIONS 

The aim of the present work is the measurement of 
the ^'^°Ba(Q;, 7)^'^'*Ce reaction cross section. The relevant 
energy range (the Gamow window) for this reaction at a 
typical 7 process temperature of 2 GK is between 5.3 and 
8.1 MeV [13] ■ In this energy range the {a, 7) cross section 
is very low, the calculations with the NON-SMOKER 
code yield values between 10~^^ and 10~^ barns. This is 
unfortunately not measurable in the laboratory, thus the 
cross section must be measured at higher energies. De- 
termined by the cross section values, the background con- 
ditions, and other experimental aspects (see below), the 
cross sections have been measured at 11 energies in the 
center-of-mass energy range between 11.6 and 16 MeV. 

Alpha capture on ^"^"Ba leads to ^'^''^Ce which is ra- 
dioactive and decays to ^"^^La with a half-life of 75.9 h. 
Subsequently, ^'^^La decays to ^■^''Ba with a half-life of 
6.67 m. This sequence allows the cross section to be mea- 
sured with the activation technique which already has 
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been used successfully for measurements relevant for the 
7 process [11]. The decay of "4Ce is followed by the 
emission of 7 rays with a very low relative intensity of 
well below 1 %. On the other hand, 7 radiation of 5 % 
relative intensity accompanies the decay of ^■^^La. The 
detection of this 7 radiation has been used to determine 
the cross section. 

Besides radiative capture, several other reaction chan- 
nels are open on ^^"Ba in the studied energy range. 
Only the (Q!,n) and (a,p) reactions have, however, non- 
negligible (above 1 /ib) cross sections. The cross section 
of ^^'^Ba{a,n)^^^Ce can again be measured with activa- 
tion since both the ground state and the isomeric state of 
^^•^Ce decay with a suitable half-life, and the decays are 
followed by high intensity 7 radiation. This measurement 
has also been carried out and owing to the different de- 
cay patterns of the ground and isomeric states of ^'^^Ce, 
partial cross section leading to these states could be de- 
termined separately. 

The (a,p) reaction can in principle also be measured 
with activation, but the product isotope, ^'^^La, has a 
continuous feeding from the decay of the (a,n) reaction 
product, ^•^■^Ce. Since the (a,n) cross section is typically 
two orders of magnitude higher than the (q;,p) one, the 
production of ^'^'^La is entirely dominated by the decay 
of ^^^Ce and therefore the {a,p) cross section could not 
be separated. 

The enriched ^'^°Ba targets (see Sec. IIII Ap contained 
substantial amount of the ^"^^Ba isotope as well. This 
isotope is also listed in Ref. [l^ in the second priority 
group of the isotopes on which a capture cross section 
measurements are needed for the 7 process. Radiative 
a capture on ^'^^Ba leads to the stable ^'^^La, hence its 
cross section cannot be measured with activation. The 
(Q!,n) reaction, however, leads to the radioactive ^■^^Ce. 
The cross section measurement of the ^■^^Ba(a,n)^'^^Ce 
reaction has therefore been also carried out in the present 
work. 

The studied reactions and the decay of the reaction 
products are shown schematically in Fig. [1] The rele- 
vant decay parameters of the produced isotopes used for 
the analysis are summarized in Table ID All decay pa- 
rameters are taken from Refs. [18rt20i] with the exception 
of the half-hfe of ^^^Ce"*. The literature value of the 
i33Qgm iiaif_iifg ]^as aj]^ unusually high uncertainty and 
the analysis of our data indicated that the half-life is 
strongly underestimated. We have therefore carried out 
a dedicated half-life measurement of this isomeric state 
resulting in a much more precise value. Details and the 
result of this work has already been published elsewhere 



TABLE I. Some decay parameters of the isotopes produced 
by the studied reactions. Only those 7 radiations are listed 
which have been used for the analysis. The data are taken 
from 18 21]. 



Reaction 


Produced 
isotope 


Half-life 


E^/keV 


relative 
7 intensity [%] 


^^'^Bafa t) 


^^^Ce 


75.9 h 


only weak gammas 


"0Ba(a,7) 


i^'-'La 


6.67 min 


604.7 


5.04 ±0.20 


"«Ba(Q,n) 


133 Qg9 


97min 


79.6 


15.9 ±2.3 








97.3 


45.5 ±6.9 








557.3 


11.4±2.3 


"''Ba(a,n) 


133 Qgm 


5.326 h 


58.4 


19.3 ±0.4 








130.8 


18.0 ±0.4 








346.4 


4.17 ±0.08 








477.2 


39.3 ±0.3 








689.5 


4.13±0.12 








784.5 


9.67 ±0.25 


"2Ba(Q,n) 


i35Ce 


17.7h 


265.6 


41.8 ±1.4 








300.1 


23.5 ±0.5 



III. EXPERIMENTAL PROCEDURE 

A. Target properties 

As it is typical for the p isotopes, the natural abun- 
dance of ^•^'^Ba is very low, only 0.1%. The use of en- 
riched targets was therefore necessary. The ^^°Ba con- 
tent of the enriched material was 11.8 ± 0.2 % and it con- 
tained also 0.7 ± 0.1 % ^'^^Ba. The enriched material was 
available in BaC03 form. Targets have been prepared 
by evaporating BaCOs onto thin Al foils. Since the 
stoichiometry of the compound may change during the 
evaporation process and the reliable determination of the 
number of ^■^"Ba target atoms (i.e., the target thickness) 
is very important for the cross section measurements, 
three different methods have been used to measure the 
target thickness. The weight of the Al foils has been mea- 
sured before and after the irradiation, the energy loss of a 
particles penetrating the targets has been measured and 
the target thickness has been also measured by Ruther- 
ford Backscattering Spectrometry. Details of the target 
thickness measurements have been published elsewhere 
[2^ . Good agreement between the three methods have 
been found and a final value of 8 % could be assigned 
to the uncertainty of the number of target atoms. The 
number density of the ^'^"Ba atoms was between 3 • 10^^ 
and 2 • 10^^ atoms/cm^. 



B. Irradiations 

The irradiations have been carried out at the cyclotron 
accelerator of ATOMKI. The targets have been put into 
the target chamber which had been used in several acti- 
vation experiments before (see e.g. [23jV 

The maximum beam current was restricted to about 
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FIG. 2. Typical 7 spectrum taken on a target irradiated by a 
14 MeV a beam shortly after the end of the irradiation. The 
peaks used for the cross section determination are indicated 
by arrows. 



1 /xA in order to avoid target deterioration. Nevertheless, 
the target stability has been continuously monitored by 
detecting the backscattered alphas with a particle detec- 
tor built into the chamber. 

The length of the irradiations varied between several 
hours and two days. Longer irradiations have been used 
at the lowest energies where the cross sections are low 
in order to increase the number of long lived ^^''Ce nu- 
clei. The half-life of the -'^■^°Ba(Q!,n) reaction products are, 
however, short compared to the length of the irradiation. 
Therefore, the time dependence of the beam current dur- 
ing the irradiation must be known. For this purpose the 
beam current was integrated and the number of current 
integrator counts was recorded in multichannel scaling 
with a time basis of one minute. 




590 600 610 620 630 640 650 
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FIG. 3. (color online) The relevant part of the 7 spectrum 
measured on a target irradiated by a 14 MeV a beam. The red 
area and label indicate the peak from the decay of ^^"^La used 
for the ^^°Ba(a,7)^''*Ce cross section measurement. Some 
other peaks from radioisotopes produced on target impurities 
are also indicated. 
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FIG. 4. (color online) The measured and fitted absolute ef- 
ficiency of the HPGe detector. In the 1 cm geometry ^Be, 
^^Mn, ''^Zn, and ^'^^Cs single line sources have been used in 
order to avoid true coincidence summing. In the 10 cm geom- 
etry ^^Na, '^■'Co, ^°Co, "^Ba, ^^^Eu, and ^^^Am sources have 
been used additionally. 



C. Detection of the induced 7 radiation 

After the end of the irradiations the targets have been 
transported to the low-background counting setup of 
ATOMKI. This setup consists of a 100% relative effi- 
ciency Canberra HPGe detector and a complete Att low 
background shielding 24?] . 

Owing to the largely different half-lives of the reaction 
products of the different investigated reactions, two sepa- 
rate 7 countings have been carried out. The first counting 



started typically 30 minutes after the end of the irradia- 
tion. In this phase the decay of the ^^^Ce^, ^^^Ce™, and 
^■^'^Ce isotopes could be measured. 

Figure [2] shows a typical 7 spectrum taken in this first 
counting period on a target irradiated by a 14 MeV a 
beam. The 7 lines used for the analysis, as listed in Ta- 
ble [H are indicated by arrows. Many other peaks are 
visible in the spectrum, they belong either to the decay 
of the studied isotopes but not used in the analysis, to 
other isotopes induced by reactions on heavier Ba iso- 
topes, or to radioisotopes produced on target impurities. 
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The spectra were stored regularly, every 10 minutes in 
the first five hours and every hour afterwards. The decay 
of the different isotopes was followed in this way and a 
half-life analysis guaranteed that no short- or long-lived 
isotopes gave contribution to the studied peaks. 

After several days of cooling, the 7 spectra of the 
targets have been taken again in order to measure the 
^^°Ba(a,7)^^''Ce cross section through the decay of ^'^"'La. 
The relevant part of the spectrum measured on a target 
irradiated at 14MeV can be seen in Fig. [3] The length of 
the counting was four days and it started approximately 
six days after the end of the irradiation. The peak of 
the ^^'^La. decay used for the ^'^°Ba(a,7)^^''Ce cross sec- 
tion determination is indicated by the red area and la- 
bel. Some other peaks are also indicated in the figure. 
They correspond to the decay of ^''As and ^''■^Ru, which 
have been produced most likely by the ^^Ga(a,n)^^As 
and ^°°Mo(a,n)^'^'^Ru reactions, respectively, induced on 
target impurities. 

Three different counting geometries have been used in 
the course of the measurements. The samples have been 
put either at 1cm, 10 cm or 27 cm distance from the end 
cap of the detector in its geometric axis. The absolute 
efficiency of the detector has been measured in all three 
counting geometries with several calibrated radioactive 
sources. The used sources were the following: ^Be, ^^Na, 
54Mn, "Co, 6OC0, 65zn, issga, i37Cs, ^^^En, ^^Am. 
Owing to the large true coincidence summing effect at 
close geometries, only single line calibration sources (^Be, 
^^Mn, ^^Zn, ^3^Cs) have been used at the 1 cm geometry. 

Figure S] shows the measured efficiencies. At 10 and 
27 cm, the measured points have been fitted with a sixth 
order logarithmic polynomial. At 1 cm the measured 
point fall into the energy region where the efficiency curve 
can well be described with a log-log straight line. The 
fitted curves can also be seen in the figure. 

Since the decay of the ^•^*'Ba(a,n)^3'^Ce reaction prod- 
ucts involves a high number of mainly low-energy 7 tran- 
sitions, the first countings have been carried out in the 
10 cm geometry in order to avoid large true coincidence 
summing effects. Nevertheless, even at this distance the 
summing effect cannot be neglected and must be deter- 
mined. This was carried out by measuring the spectrum 
of a high activity ^•^■^Ce source (produced in a high-energy 
a irradiation of a ^^^Ba target) both in the 27 and 10 cm 
geometries. At 27 cm the summing effect is well below 
1 % and therefore neglected. By taking into account the 
time elapsed betwen the two countings, a summing cor- 
rection factor for all the studied transitions was calcu- 
lated. This corrections factor was always below 10%. 

In order to maximize the count rate for the low in- 
tensity 605 keV transition of the ^^*La decay, the sec- 
ond countings have been carried out in the 1 cm geom- 
etry. In principle, the decay of ^"^^La is free from true 
coincidence summing, since the 605 keV 7 line repre- 
sents the only strong transition. The nucleus -"^^^La is, 
however, a positron emitter. Therefore, the true coin- 
cidence between the 511 keV positron annihilation line 



TABLE II. 


Measured cross section of the 


""Ba(a,7)" Ce re- 


action 






Ebcam 


Ec.m. 


cross section 


MeV 


MeV 


/ibarn 


12.0 


11.613 ± 0.045 


65.5 ± 9.8 


12.5 


12.133 ± 0.052 


118.8 ± 12.7 


13.0 


12.596 ± 0.041 


328 ± 48 


13.5 


13.080 ± 0.042 


410 ± 48 


14.0 


13.559 ± 0.046 


940 ± 67 


14.5 


14.024 ± 0.060 


1453 ± 108 


15.0 


14.509 ± 0.061 


2051 ± 187 


15.5 


14.995 ± 0.061 


2051 ± 146 


16.0 


15.509 ± 0.048 


1937 ± 142 


16.5 


15.998 ± 0.049 


1567 ± 114 



and the 605 keV line cannot be neglected. A similar 
procedure has thus been followed as in the case of the 
^'^'^Ba(a,n)^'^'^Ce measurements. Several strong ^'^^La 
sources have been prepared by the ^'^^Ba(p,n)^'^'*La reac- 
tion. For this purpose, natural BaCOs and BaCl2 targets 
have been irradiated by 8 MeV protons at the cyclotron of 
ATOMKI. The 7 spectra at 10 and 1 cm have been mea- 
sured and a correction factor of 30% ±1.3% has been 
obtained. 



IV. CROSS SECTION RESULTS 

Tables Hn mil and II VI show the measured cross sections 
of the three investigated reactions. In the case of the 
^■^'^Ba(a,n)^^'^Ce reactions, partial cross sections for the 
population of the ground and isomeric states have been 
determined separately. These are listed in Table IIIII as 
well as the total cross section. 

The effective center-of-mass energies indicated in the 
second columns have been calculated taking into account 
the beam energy loss in the target. From the known tar- 
get thickness the energy loss has been calculated with 
the SRIM code [2^. Depending on the target and the 
beam energy the energy loss was between about 20 and 
150 keV. Since the cross section is a smoothly varying 
function of the energy in the studied energy range and 
the cross section does not change much within the thick- 
nesses of the target, the effective energy has been chosen 
to correspond to the middle of the target. For the uncer- 
tainty of the effective energy the quadratic sum of half 
of the target thickness and the uncertainty of the beam 
energy has been taken, the latter one being 0.3 %. 

In the case of those reactions where more than one 
7 transition have been analyzed the results were always 
in good agreement. The error weighted mean has been 
calculated and adopted as the final cross section result. 
The uncertainty of the cross section values has been cal- 
culated using error propagation based on the following 
components: counting statistics (<15%), decay param- 
eters (< 15%), number of target atoms (8%) detec- 
tor efficiency (5%), summing correction factor (<5%), 
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TABLE III. Measured cross section of the ^^°Ba(Q,n)^^^Ce reaction 







ground state 


isomeric state 


total 


Ebcam 


Ec.m. 


cross section 


cross section 


cross section 


MeV 


MeV 


/ibarn 


/ibarn 


/ibarn 


12.5 


12.051 ± 0.083 


411 ± 82 


249 ± 36 


660 ± 116 


13.0 


12.596 ± 0.041 


983 ± 145 


611 ± 68 


1594 ± 206 


13.5 


13.080 ± 0.042 


2104 ± 310 


1620 ± 181 


3724 ± 469 


14.0 


13.559 ± 0.046 


6960 ± 1013 


5519 ± 613 


12479 ± 1549 


14.5 


14.024 ± 0.060 


14502 ± 2111 


10905 ± 1212 


25408 ± 3164 


15.0 


14.509 ± 0.061 


23762 ± 3459 


22563 ± 2506 


42651 ± 5211 


15.5 


14.995 ± 0.061 


34611 ± 5041 


36151 ± 4016 


70762 ± 8584 


16.0 


15.509 ± 0.048 


37830 ± 5511 


43548 ± 4837 


81378 ± 9803 


16.5 


15.998 ± 0.049 


37941 ± 5524 


44203 ±4911 


82144 ± 9892 



TABLE IV. Measured cross section of tiie ^'■^^Ba.{a,n)^'-^^Ce 



reaction 


Ebcam 


F, 

^c.m. 


cross section 


MeV 


MeV 


^barn 


12.5 


12.051 ± 0.083 


815 ± 159 


13.0 


12.596 ± 0.041 


1800 ± 435 


13.5 


13.080 ± 0.042 


4520 ± 923 


14.0 


13.559 ± 0.046 


11900 ± 2194 


14.5 


14.024 ± 0.060 


21500 ± 3959 


15.0 


14.509 ± 0.061 


55100 ± 10115 


15.5 


14.995 ± 0.061 


82000 ± 15053 


16.0 


15.509 ± 0.048 


95400 ± 17513 


16.5 


15.998 ± 0.049 


111000 ± 20180 



beam current integration (3%). The uncertainty of the 
"'^^^Ba(a,n)^^^Ce cross section is somewhat higher ow- 
ing to the large uncertainty (14%) of the ^'^^Ba isotopic 
abundance in the enriched ^'^'^Ba material. 
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FIG. 5. Comparison of experimental cross sections for 
^^''Ba(a,7)^"^*Ce and theoretical predictions with the codes 
NON-SMOKER [H, [13, SMARAGD jj, i28i, and TALYS 
[29| (using their default settings). 



V. IMPACT ON PREDICTIONS OF 
ASTROPHYSICAL REACTION RATES 

The experimental cross sections for the reaction 
-'^^°Ba(Q;,7)"'^'^''Ce are compared to theoretical predictions 
with the default results from the Hauser-Feshbach statis- 
tical model codes NON-SMOKER [26, 27], SMARAGD 
(version 0.8.3s) [16, 28], and TALYS (version 1.2) 29] in 
Fig. [51 The NON-SMOKER and SMARAGD predictions 
are too high by factors of 2 — 3 within the shown energy 
range. Due to additional reaction mechanisms, such as 
possible pre-equilibrium particle emission, the results of 
these codes cannot be compared to the data above 14 
MeV for this reaction but this is not relevant for the as- 
trophysical conclusions, as discussed below. The TALYS 
prediction seems to be closer to the data in the upper part 
of the shown energy range but its energy dependence is 
different from the one shown by the data, leading also to 
an overprediction of almost a factor of two at the lowest 
energy. 

It is not possible to draw a straighforward conclusion 
on the prediction of astrophysical reaction rates from this 
comparison. A closer inspection of the sensitivity of the 



resulting cross sections on a variation of the calculated 
transmission coefficients (averaged widths) is necessary 
to understand the source of the discrepancies and its 
possible impact on the reaction rate. This sensitivity is 
shown in Fig. [5] for a factor two variation of the widths. 
A sensitivity of zero implies that the cross section does 
not change when varying the width, a factor of unity 
means that the cross section changes by the same fac- 
tor as the variation factor used for the width [l^ . In the 
measured energy range, the cross sections are sensitive to 
three widths, the neutron-, 7-, and a width, respectively. 

From the a capture measurement alone it would not 
be possible to further constrain the origin of the devia- 
tion from the data because of the sensitivity to several 
widths. The combination with the (a,n) data in the same 
energy range, however, allows to check the validity of the 
a width because the (a,n) cross section is mostly sen- 
sitive to this width at these energies. The comparison 
of calculated to measured cross sections for the reaction 
^^°Ba(a,n)^^'^Ce is shown in Fig. [7] As can be seen, there 
is good agreement with the SMARAGD prediction. This 
implies that the a transmission coefficient is predicted 
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FIG. 6. Sensitivity of the cross sections for "°Ba(a,7)"*Ce 
on variation of the neutron-, proton-, a-, and 7 width. 
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FIG. 7. Comparison of experimental cross sections for 
^'^''Ba(Q,n)^^^Ce and theoretical predictions with the codes 
SMARAGD [3,[2l and TALYS Q (using their default set- 
tings). 




0.1 ^ ' ' ' ' ' ' ' 1 

23456789 10 

FIG. 8. Ratios of astrophysical reaction rates for 
^^°Ba(Q,7)^^^Ce obtained with the SMARAGD and TALYS 
codes to the standard NON-SMOKER rate 27] as function 
of astrophysical temperature Tg (given in GK). 



well and is not the reason for the overpredicted (0,7) 
cross section. The code TALYS also predicts the abso- 
lute value of the {a,n) cross section well in this energy 
range although it is using a different optical a-l-nucleus 
potential. The energy dependence, however, is different 
than the one of the data and thus it seems fortuitous 
that the absolute value is reproduced within the shown 
energy range. The deviation in energy dependence, espe- 
cially towards lower energies, is less strong than the one 
found in the (a, 7) case, implying that also the energy 
dependence of the 7- and/or neutron width may not be 
correctly described at low energy. A different a width 
energy dependence will lead to quite different cross sec- 
tions, compared to the SMARAGD prediction, at lower 
- astrophysically important - energies. 

Combining the conclusions from the (a, 7) and {a,ii) 
data it is obvious that the difference between the 
SMARAGD and TALYS predictions for the a cap- 
ture stems from differently predicted neutron- or/and 7 
widths. With the given data it is not possible to further 
constrain the uncertainties. But are these widths actu- 
ally astrophysically relevant? The astrophysical energy 
window for i30Ba(a,7)"4Ce is about 5-8 MeV Fig. 
[6] shows that the a width is the only determining factor 
in this energy range, as it is by far smaller than the 7 
width. Therefore the uncertainty in the prediction of the 
a width will dominate the resulting uncertainty in the 
prediction of the astrophysical reaction rate. 

A comparison of the astrophysical rates obtained with 
the various codes is presented in Fig. [5] Since the major- 
ity of astrophysical simulations make use of reaction rates 
from [13 which were obtained with the NON-SMOKER 
code, the ratio of the more recent predictions to this stan- 
dard is plotted. The predictions differ by not more than 
20% at the highest temperatures but they run apart at 
lower temperatures f3(^. The relevant temperature re- 
gion is 2 — 3 GK, at the lower end of the shown range. 
It is noteworthy that these results were obtained with 
the standard settings of the codes, not attempting to fit 
any further energy dependence of the optical a-l-nucleus 
potential to the data. It has to be realized that even the 
reproduction of the energy dependence of the averaged a 
width in the measured energy range - as provided by the 
potential by [3l| used with the SMARAGD and NON- 
SMOKER codes here - does not guarantee that the en- 
ergy dependence is similarly well described at much lower 
energies. This underlines the fact that it is necessary to 
measure within the astrophysically relevant energy win- 
dow in order to constrain the reaction rates. 

A very similar picture arises when comparing the mea- 
sured -'^^^Ba(a,n)^^^Ce cross sections to predictions, as 
shown in Fig. O The theoretical calculations reproduce 
the data well (with similar deviations from the measured 
energy dependence as in the ^'^°Ba case) which leads to 
the conclusion that the a widths are described well in 
the measured energy range. Also the comparison of the 
actual stellar reaction rates for ^■^^Ba(Q;,7)^'^^Ce obtained 
with the different codes in Fig.[Tn]looks similar to the pre- 
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FIG. 9. Comparison of experimental cross sections for 
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SMARAGD [ll,[2i| and TALYS H (using their default set- 
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FIG. 10. Ratios of astrophysical reaction rates for 
^^2Ba(Q,7)^^^Ce obtained with the SMARAGD and TALYS 
codes to the standard NON-SMOKER rate [IJ as function 
of astrophysical temperature Tg (given in GK). 



vious case: the newer predictions are lower at high tem- 
perature than the NON-SMOKER rates but show a dif- 
ferent temperature dependence. While the SMARAGD 
rate is considerably higher in the relevant temperature 
range, the TALYS rate is comparable or lower than the 
standard rate. 



VI. CONCLUSIONS 

The cross sections for the reactions ^■^°Ba(Q;,7)^'^'*Ce, 
^^''Ba(a,n)^^^Ce, and ^^^Ba(Q!,n)^^^Ce have been mea- 
sured closely above the (a,n) threshold with the acti- 
vation technique. The combination of (01,7) and (a,n) 
data allowed to test the predicted a widths (and thus 
the optical a-l-nucleus potential) in the measured energy 



range. Good reproduction of the ^^"Ba(a,n)^^^Ce and 
^'^^Ba(Q!,n)^'^^Ce cross sections was found when using the 
potential by [sH , which is a standard potential often used 
in astrophysical applications. 

Assuming the validity of the potential by [3lj also at 
low energies, we present Table FVl with the stellar reactiv- 
ities obtained with the SMARAGD code, version 0.8.3s, 
for the ^'^'^Ba(a,7)^'^'*Ce reaction and its inverse. The 
parameters for a fit of the reactivities in REACLIB for- 
mat are given in Table IVII Stellar reactivities include 
transitions from thermally populated states of the target 
nuclei. Even if it were possible to measure in the as- 
trophysically relevant energy range, only the transitions 
from the ground state of the target would be determined. 
The ground state contribution X to the stellar rate as in- 
troduced in Ref. U^] is also given, showing by how much 
the uncertainty may be reduced by a measurement on 
-'^'^"Ba or ^■^*Ce, respectively, in the ground state. Com- 
paring Xa--^ for the capture reaction with X~^a for the re- 
verse (7,0;) reaction also nicely shows how much smaller 
the measured contribution to the stellar rate would be 
when using a photodisintegration experiment instead of 
a capture one. 

The reaction ^■^*Ce(7,Q;)^'^"Ba received special mention 
in [1^ because its rate is close (within a factor of 10) to 
the "'^■^'*Ce(7,n)^'^'^Ce rate at plasma temperatures from 
2 — 2.5 GK. The photodisintegration path depends on 
the competition between charged particle and neutron 
emission. Due to the large uncertainty in the prediction 
of the low-energy a width, the path may not be well de- 
termined when the (7,0?) and (7,n) rates are within such 
a factor. At 2 GK, the (7,0) rate was almost a factor 
of 10 faster than the (7,n) rate in [l^, at 2.5 GK it was 
already lower by about 1/5. With the new rates from Ta- 
bleEl the (7,a)/(7,n) rate ratio of "4Ce is 100 at 2 GK 
whereas the two rates are comparable at 2.5 GK. The dif- 
ference comes from the newly calculated ^'^*Ce(7,a)"'^^°Ba 
rate because the SMARAGD and NON-SMOKER pre- 
dictions of the ^■^''Ce(7,n)^'^'^Ce rate are almost the same. 
The low-energy difference in the (7, a) rates comes from 
an improved numerical treatment of the charged-particle 
transmission coefficients at low energy in SMARAGD 
(see section 5.4.2 in [l6j| for a discussion). Thus, the 
path deflection at 2 GK is not shifted with respect to the 
one found in (Tsj but a branching is introduced at 2.5 
GK which removes some flux from the (7,n) path lead- 
ing to proton-richer Ce isotopes. Previously, the flrst 
path deflection point (where (7,0;) dominated over (7,n)) 
at 2.5 GK was at ^'^^Ce. It may lose importance with 
the new rates because part of the photodisintegrations 
branch off already at ^'^'*Ce. This implies that the 7 pro- 
cess path would lie slightly closer to stability. It has to 
be cautioned, however, that the rates in the astrophysical 
energy window are still not well constrained by experi- 
ment. 

As mentioned in Sec. |V] the deviation of the predictions 
from the data in Fig. [5] is not caused by a mispredicted a 
width but by problems in the predicted neutron and/or 
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TABLE V. Stellar reactivities R in cm^s ^mole ^ and ground state contributions X to the stellar rate as function of temperature 
T. Both quantities are given for the reactions ^'^'^'B&{a,'yy^'^ Ce and Ce{'y ^aY^^Iiau. 



T (GK) 


















0.15 


1.173 


X 


io-»* 


1.0 


5.4699 


X 


10-''^ 


0, 


.1882 


0.20 


1.137 


X 




1.0 


8.0374 


X 


10-63 


0, 


.1787 


0.30 


1.363 


X 


10-58 


1.0 


1.7428 


X 


^0-48 


0, 


.1501 


0.40 


7.683 


X 


10-5° 


0.9998 


1.5008 


X 


10-39 


0, 


.1318 


0.50 


1.094 


X 


io-« 


0.9988 


2.9757 


X 


10-33 


0, 


.1197 


0.60 


4.804 


X 


10-39 


0.9950 


1.7158 


X 


10-2« 


0, 


.1105 


0.70 


2.235 


X 


10-35 


0.9868 


1.0086 


X 


10-24 


0, 


.1032 


0.80 


2.204 


X 


10-32 


0.9727 


1.2250 


X 


10-21 


0, 


,09701 


0.90 


7.044 


X 


10-30 


0.9522 


4.7184 


X 


10-19 


0, 


.09150 


1.00 


9.621 


X 


^0-28 


0.9259 


7.6701 


X 


10-" 


0, 


.08637 


1.50 


2.119 


X 


10-20 


0.7439 


3.4747 


X 


10-9 


0, 


05929 


2.00 


5.827 


X 


io-« 


0.5558 


1.6807 


X 


10--* 


0, 


.03129 


2.50 


6.798 


X 


10-13 


0.4032 


3.0934 


X 


10-1 


0, 


.01398 


3.00 


1.378 


X 


10-10 


0.2911 


9.1856 


X 


lOi 


n 




3.50 


9.264 


X 


10-9 


0.2126 


8.6161 


X 


103 


0, 


.002428 


4.00 


2.783 


X 


10-^ 


0.1594 


3.4934 


X 


105 


0, 


.001085 


4.50 


4.397 


X 


10-6 


0.1230 


7.1454 


X 


10® 


5, 


.396 X 10-* 


5.00 


4.096 


X 


10-5 


0.09718 


7.8862 


X 


10^ 


2, 


.971 X lO-"* 


6.00 


1.010 


X 


10-3 


0.06354 


1.8376 


X 


109 


1, 


.096 X 10-* 


7.00 


7.134 


X 


10-3 


0.04276 


1.0434 


X 


lOlo 


4, 


.365 X 10-5 


8.00 


2.278 


X 


10-2 


0.02889 


3.1260 


X 


lOlo 


1, 


,573 X 10-5 


9.00 


4.550 


X 


10-2 


0.01923 


6.5428 


X 


lOlo 


4, 


,662 X 10-6 


10.00 


6.855 


X 


10-2 


0.01250 


1.0774 


X 


10" 


1, 


,168 X 10-6 



TABLE VI. Fit parameters for the REACLIB fit format [27'| for i30Ba(a,7)i3*Ce and its inverse. Note that the value for 
i3*Ce(7,a)i30Ba obtained in this parameterization has to be multiplied by the ratio of the temperature-dependent partition 
functions G(i30Ba)/G(i3'*Ce) 2713 to get the stellar reactivity. 





ao 


ai 


a2 


as 


04 


as 


as 


(a,7) 


-842.1489 


-75.13541 


-1326.36 


2286.821 


-112.6802 


5.078861 


-1117.492 


(7,a) 


-817.10229989 


-75.12610319 


-1326.36 


2286.821 


-112.6802 


5.078861 


-1115.992 



7 width. Although this is inconsequential regarding the 
astrophysical -'^^''^Ce(7,a)^^°Ba photodisintegration reac- 
tion, it may affect the prediction of the (7,a)/(7,n) rate 
ratio discussed above. The ^^'^Ce(n,7)^^''Ce rate (and its 
inverse ^•^*'Ce(7,n)^'^^Ce) mainly depends on the neutron 
width at 2 GK and is sensitive to both neutron and 7 
width at temperatures 2.5 and 3 GK. Unfortunately, the 
dependences in the measured ^'^*'Ba(Q;,7)^'^*'Ce cross sec- 
tion and the inferred ^^'^Ce(n,7)^'^*Ce cross section are 
different. The overpredicted ^■^*'Ba(a,7)^'^*'Ce cross sec- 
tion can be due to a 7 width being too large or a neutron 
width being too small. On the other hand, a decreased 
7 width would lead to a decrease in the predicted rate 
for i33Ce(n,7)"4Ce at T < 2.5 GK while an increased 
neutron width leads to an increase in the same rate for 
T < 3 GK. Without further information regarding the 
neutron and 7 widths in ^^"^Ce it is not possible to draw 
further conclusions. 

Similarly to the ^^°Ba(Q;,7)^^'*Ce case discussed 
above, we also present the stellar reactivities for 
^^^Ba(Q;,7)^^^Ce and its reverse rate in Table IVlIl Again, 



the potential of 31] is used with the SMARAGD code, 
version 0.8.3s. The parameters for a fit of the reactivi- 
ties in REACLIB format are given in Table IVllIl along 
with the ground state contributions Xa-y and X^a to the 
stellar rate for the (a, 7) and the (7, a) rate, respectively. 
The X^fo. for this reaction are comparable to the ones for 
i34Ce(7,a)"0Ba. 

The reaction ^'^^Ce(7,Q!)^'^^Ba was pointed out in [l^ 
because ^^^Ce could become a new deflection point in 
the photodisintegration path if this rate is increased. 
As explained above, this rate has to be compared to 
the ^■^^Ce(7,n)t35Ce rate which is also only theoretically 
known. Again, the prediction of the (7,n) rate remains 
almost unchanged with respect to the NON-SMOKER 
value. Therefore, the much higher new (7,0) rate ex- 
ceeds the (7,n) rate already at low plasma temperature 
and therefore the (7,n) flow from stability would be in- 
terrupted already at this isotope. This renders the above 
discussed ^'^^Ce deflection (or branching) unimportant 
and forces the 7 process flow to move even closer to stabil- 
ity. It implies that the p isotope t-^OBa does not directly 
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TABLE VII. Stellar reactivities R in cm^s ^mole ^ and ground state contributions X to the stellar rate as function of tem- 
perature T. Both quantities arc given for the reactions ^^^'Ba.{a,'yy'^^ Ce and ^'^^Ce('y,ay'^^'B&. 



T (GK) 


















0.15 


1.156 


X 




1.00 


2.0268 


X 


10-90 


0, 


.1822 


0.20 


1.085 


X 




1.00 


2.0713 


X 


^0-74 


0, 


.1709 


0.30 


1.303 


X 




1.00 


3.2317 


X 




0, 


.1425 


0.40 


7.324 


X 


10-50 


1.00 


2.3519 


X 


^Q-45 


0, 


.1254 


0.50 


1.046 


X 


io-« 


0.9999 


6.7055 


X 


10-38 


0, 


.1136 


0.60 


4.611 


X 


10-39 


0.9994 


2.2899 


X 


10-32 


0, 


.1042 


0.70 


2.155 


X 




0.9977 


4.7843 


X 


10-28 


0, 


.09549 


0.80 


2.138 


X 


10-32 


0.9941 


1.5035 


X 


10-24 


0, 


,08691 


0.90 


6.882 


X 


10-30 


0.9876 


1.2147 


X 


10-21 


0, 


.07838 


1.00 


9.487 


X 




0.9777 


3.5622 


X 


10-19 


0, 


,07006 


1.50 


2.248 


X 


10-20 


0.8766 


9.5831 


X 


10-11 


0, 


,03618 


2.00 


6.715 


X 




0.7316 


1.1248 


X 


lo-'^ 


0, 


,01739 


2.50 


8.242 


X 


10-13 


0.5934 


3.4109 


X 


10-2 


0, 


,008395 


3.00 


1.642 


X 


10-10 


0.4821 


1.3169 


X 


lOi 


0, 


,004211 


3.50 


9.888 


X 


10-9 


0.3993 


1.3340 


X 


103 


0, 


,002275 


4.00 


2.457 


X 


10-^ 


0.3379 


5.0494 


X 


10-1 


0, 


,001349 


4.50 


3.043 


X 


10-® 


0.2895 


8.6088 


X 


lO'^ 


8. 


,633 X 10-* 


5.00 


2.151 


X 


lo-'^ 


0.2482 


7.4991 


X 


10® 


5, 


,756 X 10-1 


6.00 


3.090 


X 


10--* 


0.1781 


1.3316 


X 


108 


2, 


,509 X 10-* 


7.00 


1.534 


X 


10-3 


0.1208 


7.9416 


X 


108 


8, 


,864 X 10-5 


8.00 


4.270 


X 


10-3 


0.07657 


2.7269 


X 


109 


2, 


,342 X 10-5 


9.00 


8.507 


X 


10-3 


0.04542 


6.6817 


X 


109 


5. 


,033 X 10-® 


10.00 


1.358 


X 


10-2 


0.02575 


1.2869 


X 


lOlo 


9, 


,795 X 10-^ 



TABLE VIII. Fit parameters for the REACLIB fit format ^ for i32Ba(a,7)i3®Ce and its inverse. Note that the value for 
i3®Ce(7,a)i32Ba obtained in this parameterization has to be multiplied by the ratio of the temperature-dependent partition 
functions G(i32Ba)/G(i3®Ce) ^ to get the stellar reactivity. 





ao 


ai 


a2 


as 


04 


as 


as 


(0,7) 


-1155.448 


-121.1251 


25.24542 


1271.057 


-87.73028 


4.764612 


-404.0757 


(7, a) 


-1130.40071984 


-126.44344235 


25.24542 


1271.057 


-87.73028 


4.764612 


-402.5757 



receive contributions from the Ce isotopic chain but only 
via (7,n) reaction chains from neutron-richer Ba isotopes. 
The same cautionary remark as made above, however, 
also applies here. The validity of the potential by [3l| 
has been proven by the experiment only at higher en- 
ergy than required in the calculation of the astrophysical 
reaction rate. 

In summary, the present experiment and its interpre- 
tation are exemplary for the typical difficulties encoun- 
tered in the determination of charged-particle cross sec- 
tions and reaction rates for astrophysics and, more specif- 
ically, the p process. The cross sections at astrophysical 
energies are unmeasurably tiny. On the other hand, the 
cross sections at higher, experimentally accessible, ener- 
gies show different sensitivities to nuclear properties than 
those at lower energy because more reaction channels are 
open and more particle widths contribute. A further fun- 
damental difference is that charged-particle widths at as- 
trophysically relevant energies tend to be much smaller 
than the 7 width, which is not the case at the usually 
investigated higher energies. Regarding a capture, the 



leading uncertainty in the astrophysical rate comes from 
the optical a+nucleus potential. Measuring at higher 
energies, though, the various contributions and depen- 
dences of the particle widths and the 7 width have to be 
disentangled by combining data for different reactions. 
Regardin g th e a width, this can be achieved by (n,a) 
reactions [32|, but also by (a,n) measurements which un- 
fortunately have a strongly negative reaction Q value in 
proton-rich nuclei. Therefore the optical potential can 
only be tested at comparatively high energies and the- 
oretical considerations have to be invoked to arrive at 
an astrophysical rate. Nevertheless, {a,n) measurements 
even close to the neutron emission threshold are scarce 
for this mass region. Further (a,n) data are required and 
also further experimental efforts to access (a,7) cross sec- 
tions at even lower a energy. 
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